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Abstract 

We study the influence of two-loop radiative corrections of en- 
hanced electroweak strength on Higgs production at the LHC. We 
consider Higgs production by the gluon fusion mechanism, with the 
subsequent decay of the Higgs boson into a pair of Z bosons, and in- 
corporate the resonance shape corrections up to order {g^ m^j / ni^)"^ ■ 
We take into account the full gg — > ZZ process and the qq ZZ 
background, as well as the subsequent decay of the Z pair into leptons. 
We also discuss the theoretical uncertainty related to the use of the 
equivalence theorem in this process. 



1 Introduction 



The arguably most important task of the Large Hadron Colhder is to pro- 
vide evidence regarding the nature of the mechanism responsible for the 
spontaneous breaking of the electroweak symmetry. This is crucial for un- 
derstanding the physics at the electroweak scale and may give hints about 
the dynamics underlying the standard model at higher energies. 

The simplest model of spontaneous symmetry breaking, the Higgs sector 
of the minimal standard model, will be testable at the LHC up to an energy 
scale of ~ 1 TeV. This covers the mass range where the Higgs particle is 
believed to make sense as a fundamental field. This is related to the triviality 
of the Higgs sector, which is a nonperturbative problem. Lattice simulations 
of the (j)^ theory seem to indicate an upper bound of the order of 650 GeV 
for rriH, beyond which the Higgs mass would be larger than the cutoff scale 
[^. Therefore a strongly interacting Higgs sector would not exist. However, 
these results must be interpreted with caution. There are indications that 
they depend on the lattice regularization, and in fact the triviality bound 
may be considerably higher [^. Also the presence of a large Yukawa coupling 
may play a role because it relaxes the bound derived from the position of the 
Landau pole. Unfortunately, it is difficult to check this nonperturbatively 
because of notorious problems with treating fermions on the lattice in the 
presence of Yukawa interactions Q. One cannot exclude the possibility of a 
Higgs mass of the order of 1 TeV. 

The Higgs boson is produced in hadron collisions via the gluon fusion 
and the vector boson fusion processes. With mt ~ 180 GeV, the gluon fusion 
process dominates unless the Higgs is very massive, of the order of 1 TeV [Q]. 
Heavy Higgs bosons mainly decay into pairs of vector bosons. The decay 
channel H ZZ — > 2{l^l~) provides a clean signal if mn is smaller than 
~ 800 GeV. For a heavier Higgs, the channel H — > ZZ — > l'^l~vv must be 
included. The Higgs then appears as a Jacobian peak in the pT distribution 
of the Z bosons. To approach the 1 TeV mass range, one needs to consider 
the H — > ZZ — > + 2 jets and H W^W~ channels, which are affected 
by large backgrounds. 

The Higgs production mechanisms at hadron colliders were studied ex- 
tensively at leading order Q. The NLO QCD corrections to the gluon fusion 
process were calculated recently, and indicate a substantial increase of the 
Higgs production |5|, ^ . It is the purpose of this paper to study the NNLO 
corrections of enhanced electroweak strength to Higgs production via gluon 
fusion, with the Higgs subsequently decaying into a Z pair. 
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The large mn range is affected by potentially large electroweak radiative 
corrections due to the strong selfinteraction of the Higgs sector. For instance, 
the radiative corrections of increased electroweak strength to the Higgs decay 
into vector bosons are known at two-loop order and are of the order of 
24% for rriH ~ 900 GeV. It is therefore legitimate to ask what is the effect 
of such corrections on the production and decay of heavy Higgs bosons. 

A consequence of the Slavnov-Taylor identities of the theory is the equiv- 
alence theorem. This theorem relates in the high energy limit the amplitudes 
with external longitudinal vector bosons Vl and the corresponding ampli- 
tudes with the Vl replaced by Goldstone bosons ||8|-|11]. The equivalence 
theorem is a useful tool for calculating leading contributions in the Higgs 
mass in a simple way, but within this context the resonance region should be 
treated carefully. The equivalence theorem holds order by order in pertur- 
bation theory, while in the resonance region a Dyson summation is needed 
which contains incomplete contributions from higher orders. If one fails to 
carry out the calculation consistently at the given order in the coupling con- 
stant, the result will differ from the corresponding leading rrin contributions 
calculated in the full electroweak theory (see for instance ref. and |14], 
and references therein). 

Related to the large coupling domain is also the problem of unitarity vio- 
lation when the Higgs mass exceeds certain limits. Unitarity breaking effects 
are always present within perturbation theory and are formally of higher or- 
der in the coupling constant. Numerically, they may become important if 
the coupling constant is large, and raise questions about the relevance of 
calculations of physical processes. 

A few approaches were proposed in the literature to deal with such prob- 
lems. To mention only two of them, ref. [13| proposes the use of a momen- 
tum dependent width of the Higgs in the resonant propagator instead of 
a constant width in order to reduce the size of the unitarity violations in 
longitudinal vector boson scattering. Ref. ||l^ notes that this procedure re- 
produces the correct result in the resonance region but has a bad behavior at 
higher energies, where the nonresonant propagator should be used. There- 
fore it proposes a summation of the full tree level Goldstone boson scattering 
amplitude rather than a Dyson summation of only the Higgs selfenergy in 
order to interpolate between the two regions. 

These approaches amount in fact to changes in the incomplete contri- 
butions of higher order in the coupling constant. One should see them as 
leading order calculations. They are not complete at one-loop order, and 
their use is largely motivated by the fact that they lead to a solution which 
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has the desired numerical behaviour in the specific process considered. For 
instance, using the momentum dependent width in the Higgs resonant prop- 
agator gives a numerically well behaved amplitude in ref. |jl^, but would 
lead to an unwanted result in the case of fermion-fermion scattering. In- 
deed, the one-loop momentum dependent selfenergy used in ref. [13| would 
result in a shift in the position of the Higgs resonance in fermion scattering 
at leading order in mn, which in fact is absent in the full one-loop radia- 
tive corrections |16|. In this case the constant width is obviously the better 
choice. 

A systematic treatment would instead require a complete calculation 
of radiative corrections to the given process. This ensures both that the 
equivalence theorem is valid up to the desired order, and that the unitarity 
violating contributions are of a higher order in the coupling constant com- 
pared to the tree amplitude, which is the best one can do within perturbation 
theory. 

This paper is devoted to the study of the two-loop radiative corrections 
of enhanced electroweak strength to the gg — > ZZ process. By using the 
equivalence theorem, we derive in the following section the two-loop correc- 
tions to the shape of the Higgs resonance in the related process tt — > if — > zz. 
We then study the influence of these corrections on Higgs searches by means 
of a Monte-Carlo simulation of the full gg — > ZZ — > 4/ process and of the 
qq ZZ — > 4/ background at LHC. 



2 Resonance shape corrections 

In this section we derive the NNLO radiative corrections to the scattering 
process ti ^ H ^ zz, where z are Goldstone bosons. We are interested only 
in radiative corrections of enhanced electroweak strength, that is, corrections 
which grow like mjj at one-loop order and like rrij^ at two-loop. 

In the following, we denote by A the selfcoupling of the Higgs field, 

with g^ = 4^/2 Gp, Gp = 1.16637- IQ-^ GeV~'^, and mw = 80.22 GeV . 

The radiative corrections of leading order in niH to the ti ^ H ^ zz 
process are depicted in fig. 1. They consist of the radiative correction 
to the Yukawa coupling Vhu^ which is momentum independent ||l^, p^ , 
the momentum dependent correction to the Hzz vertex Vhzz[s), and the 
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correction to the Higgs propagator n(s). In the resonance region one needs 
to perform an all-order Dyson summation of the proper selfenergy. Far 
from the peak, this resonant propagator will differ from the nonresummed 
two-loop selfenergy insertion by contributions of order A^. 
At leading order, the diagram of fig. 1 reads simply 

(2) 



s — mjj + irriH^H 

up to the tree level Htt and Hzz couplings, which are constants independent 
of s. Fh in eq. 2 is the tree level Higgs width. This expression is correct at 
leading order, i.e., 0{X~^) on the resonance and 0{1) far from the resonance. 
The radiative corrections change this expression into: 

Vmt T^. t7t-^Vh.^{s) . (3) 

We are primarily interested in radiative corrections in the resonance 
region. We therefore perform a momentum expansion around s = m|^: 

n(.) = n(ml,) + ^^n' + i^^^n" + . . . 

vhUs) = vh.mI) + '-^yL. + 'j'^'^' yL. + ■ • • ,(4) 

H H 

For treating the radiative corrections consistently as an expansion in the 
coupling constant A, one notes that in the resonance region the quantity 
s — rn\j is of the order of mH^H-, that is, 0{X). This must be taken into 
account together with the loop expansion of the quantities V^tf, 11 (m^), 

VHzzi'm'fj), n', V^2^, To have the complete NNLO amplitude, one 

needs to keep the following contributions: 



3\ 
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The correction to the Yukawa coupHng reads [18| |19]: 



(l—loop) 



V, 



Htt 

(2-loop) 



Htt 



, 13 vr > 

( ^ A = .132325 A 

^16 8 ^ 

(-.26387 ±1.3- 10-4) A^ . 



The correction to the Hzz vertex was calculated in ref. 
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(.17505951 - a.4190989) A 

- (.53673 ± 4.1 • 10~4) + i (.32811 ± 3.1 • lO^^; 



A^) 



By evaluating the one-loop diagrams which contribute to the H ^ zz 
decay, on finds also: 



/ (l-loop) 



Hzz 



1 + 



\/3 vr 5 vr^ 



12 
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log 2 



48 

i. 15169744) A 
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Let us now turn to the correction to the Higgs propagator. 

At the order in which we are working, the quantity n(s) is the imagi- 
nary part of the proper Higgs selfenergy. The real part of the Higgs proper 
selfenergy at s = m\ and its first derivative with respect to s are absorbed 
into the Higgs mass counterterm and wave function renormalization, respec- 
tively. The real part of the selfenergy may contribute an imaginary piece to 
eq. 5 only at ©(A^), through the quantity n"(^~'°°P). Therefore n(m|^) is 
the total width of the Higgs including the appropriate radiative corrections. 

At tree level, the main decay channels of a heavy Higgs are given by: 
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The radiative corrections to the decay widths of eqns 
two-loop level in the heavy Higgs approximation [O, 19, ^: 



(9) 

9 are known at 
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+ .97103 ± 8.2 • 10 
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H^W+W- ,ZOZO 



24 4 

1 + .350119 A + f .97103 ± 8.2 • 10~^ ) A^ 
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H^tt 



-p(tree) 



-,(4ree) 
H^tt 



1 + .264650 A - f .51023 ± 2.5 • 10~^ 1 A^ 



(10) 



Note that in eqns. 10 some incomplete subleading contributions are 
present in the radiative corrections, as discussed in ref. Q. They appear 
if one multiplies the full tree level width by the radiative correction factor 
calculated in the leading ran approximation. These terms are of the same 
order in the coupling constant as the theoretical uncertainty related to the 
use of the equivalence theorem while calculating radiative corrections. It is 
thus not possible to decide unambiguously whether it is better to keep them 
or to drop them without calculating the complete subleading contributions 
explicitly. Numerically, this ambiguity is at 1% level at most. As such it 
can be safely neglected. 

With these results, one can identify in eq. 5: 



n(i-'°°p)(m^) + n(2-'°°p)(m^) + n(3-'°°p)(m^) = 

T + r^^^o^o + Tu_^ii + T H^Aw ,wwzz ,iz 1 



(11) 



with Fj^^^y+iy-, TH^zoz<i and T^^^t given by eqns. 10. 
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The term Th^4w,wwzz,4.z originates from the three-loop cut diagrams 
shown in fig. 2. All other cut diagrams of three-loop selfenergy diagrams 
are allready contained in ^^^^^^+yy- z^zo- 

Including or dropping Tu_^ii in eq. 11 is in principle irrelevant at the 
order in which we are working because this is a contribution of higher order 
in the top quark mass. In practice this is not a large effect because the ti 
branching ratio is small and because of the partial cancellation between the 
one-loop and the two-loop corrections to this partial width. 

The quantity n'(^~'°°^\ which also contributes to eq. 5, was calculated 
inref. [|l|: 



4 V 12 48 

2 



= 1.002245142 m// . (12) 
To evaluate Th^aw wwzz Az, let us introduce the following notations: 



^H-^Az = fAz A To , 

3 TT 

^H^Aw = ^Aw A To , To = -3- mnX ■ (13) 

o 

We calculate the dimensionless factors ip by means of a Monte-Carlo 
integration over the four-body phase space of the cut diagrams of fig. 2. 

The diagrams themselves were calculated through the equivalence the- 
orem, therefore only the leading mn terms are kept at the level of the 
amplitude. However, when integrating over the phase space, we allow for 
a finite mass of the vector bosons by keeping the momenta of the outgoing 
Goldstone bosons at an invariant mass equal to the vector boson masses, 
mz and mw- This procedure is similar to the way we multiplied in eqns. 
10 the full tree level width by the loop corrections obtained via the equiv- 
alence theorem. In both cases the result contains incomplete subleading 
terms from the phase space integration (eq. 10 contains also contributions 
from the transversal degrees of freedom of the vector bosons) . This is moti- 
vated by the assumption that the phase space factor is numerically the main 
subleading effect. To check the reliability of the equivalence theorem calcu- 
lation of this process, we calculated the decay H 4Z^ in the electroweak 
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fAw 


,AET) 






"PAz 


600 


4.14 • 10"^ 


1.88 • 10^^ 


8.97 • lO"'' 


9.07 • 10"^ 


1.69 • 10"^ 


700 


6.42 • 10 ^ 


3.34 • 10 ° 


1.58 • 10"^ 


1.67 • 10"^ 


4.67 • 10 ° 


oUU 


0.41/ • ±U 




00 1 n— 5 


"^s 1 n— 5 


8 7/1 1 n— 1> 


900 


1.02 • lO-'^ 


6.24 • 10-^ 


2.79 • 10-^ 


3.01 • 10"^ 


1.33 • 10"^ 


1000 


1.18 • 10""^ 


7.53 • 10"^ 


3.25 • 10"^ 


3.51 • 10"^ 


1.78 • 10"^ 


1100 


1.31 • 10"* 


8.70 • 10"^ 


3.66 • 10"^ 


3.93 • 10"^ 


2.22 • 10'^ 


1200 


1.41 • 10-* 


9.65 • 10-^ 


4.01 • 10-^ 


4.29 ■ 10-^ 


2.62 ■ 10"^ 


1500 


1.64 • 10-* 


1.20 • 10"* 


4.73 • 10"^ 


5.00 • 10-^ 


3.59 • 10"^ 


2000 


1.87- 10"* 


1.45 • 10"* 


5.41 • 10"^ 


5.61 • 10"^ 


4.61 • lO"*' 


3000 


2.05 • 10"* 


1.66 • 10"* 


5.93 • 10-^ 


6.04 • 10"^ 


5.51 • 10"^ 


5000 


2.16 • 10"* 


1.80 • 10"* 


6.24 • 10"^ 


6.29- 10~^ 


6.08 • 10'^ 


oo 


2.23 • 10-* 


1.89 • 10-* 


6.41 • 10-*^ 


6.41 ■ 10-^ 


6.41 • 10"^ 



Table 1: The value of the factors (p calculated via the equivalence theorem 
((^(^^)) and in the electroweak theory, as a function of the Higgs mass. For 
(p{ET) ^ momenta of the outgoing Goldstone bosons have invariant masses 

mz and mw- fAz'^ corresponds to the decay H — calculated in the 

electroweak theory in the approximation (livpi) ^ Vpul'^z- 'Viz '^^ exact 
result, summed over the polarization states of the Z bosons. 

theory. We calculated the decay H — ^ 4Z£ in the approximation that the 
longitudinal polarization vectors of the Z are given by €l{Ph) ~ Pfi/n^z, as 
well as the complete H AZ^ decay. 

The results arc given in table 1. The limit m^y = mz = 0, in which the 
equivalence theorem is exact, corresponds to tuh oo. One can see that 
the equivalence theorem with the phase-space factors included approximates 
well the electroweak decay width in the approximation €l{p^) p^/mz- 
However, the equivalence theorem and the exact H AZ^ result start to 
agree only for uih of the order of 2.5 — 3 TeV. The equivalence theorem is 
not a good approximation for uih ~ 1 TeV because uih is not large enough 
compared to the four-particle threshold at ~ 360 GeV. 

This is to be contrasted to the situation encountered in one- and two- 
loop calculations. At two-loop level there are only cut diagrams with two 
massles particles on the cut lines because the Goldstone bosons are always 
produced in pairs. In this case the equivalence theorem works better than 
at three-loop level, where diagrams with four massles particles on the cut 
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lines are allowed, and the threshold in the full electroweak theory is higher. 

Generally, for given masses of the vector bosons and for fixed external 
momenta, one expects the equivalence theorem to become a progressively 
bad approximation as the number of loops increases. The subleading con- 
tributions become increasingly important because of the presence of multi- 
particle cuts. 

In this calculation, the uncertainty related to subleading contributions 
in the four Goldstone decay is numerically negligible. Compared to the two 
Goldstone channel, the four Goldstone decay is strongly suppressed by phase 
space. 

We have now all ingredients needed to evaluate eq. 3 at NNLO. At 
lowest order, the expression which we obtain is equivalent to eq. 2. At the 
same time, it incorporates the full radiative corrections of 0{X) and O(A^). 
It also contains some incomplete subleading contributions. 

In fig. 3 we compare the LO expression to the NNLO calculation. 

Although the purpose of this calculation was primarily to obtain the 
radiative corrections to the shape of the resonance in a consistent way, one 
notes that the amplitude is well behaved far from the resonance as well. 

Compared to the corrections to the Higgs width, the resonance shape 
corrections of fig. 3 are relatively small. This is because of the partial com- 
pensation of the effect of the radiative corrections to n(s) and to Vhzz{s)- 
For the same reason, the shift of the resonance peak towards lower momen- 
tum is smaller than the shift in fermionic scattering |16|. 

Finally, let us notice that the NNLO is the lowest order where nontrivial 
corrections to the resonance shape occur. At NLO there is only a correction 
to the Higgs width because the momentum dependence of the quantities 11 
and Vhzz is of higher order in A. 



3 Gluon fusion at hadron colliders 

The diagrams which contribute to the process gg ZZ are shown in fig. 4. 
Among these, only the Higgs production diagram receives radiative correc- 
tions of enhanced electroweak strength. 

The leading niH radiative corrections to the Higgs production diagram 
of fig. 4 a) are identical to the corrections to the process ti ^ H ^ zz, 
which were discussed in the previous section, by the equivalence theorem. 
This is true in the limit of large top mass as well. In this limit the top 
quark does not decouple, and the gg ^ H subgraph results in the effective 
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interaction Cgjii = g Og/ {24: tt m\Y) H G^,^ G"'^'^ , which receives the same 
radiative corrections at leading order in mn as the Yukawa coupHng. 

We incorporated the radiative corrections of the previous section in a 
Monte-Carlo event generator to calculate the gg — > ZZ process and the 
qq — >■ ZZ background at LHC. The subsequent decay of the Z pair into 
fermions is included in the narrow width approximation by using the density 
matrix formalism. Details of the calculation can be found in ref. pO| , 

We take y/s = 14.5 TeV for the CM energy of the LHC, and rrit = 180 
GeV for the mass of the top quark. We use the MHRS parton distribution 
functions |^], with A = 190 MeV, evolved to the scale = s/4, and as = 
127r/(231og Q^/A^) corresponding to five flavours. As a rough simulation 
of the detector geometry, we impose a rapidity cut on the outgoing leptons 
of \yi\ < 2.5, and request that the leptons have a transverse momentum 
larger than 20 GeV. The cross sections shown in the following correspond 
to four muons in the final state. The branching ratio should be multiplied 
by a factor 4 for muons or electrons in the final state, and 24 for l^l~vD 
(with these cuts, the cross sections for neutrinos in the final state will be 
underestimated to some extent). One can consider an integrated luminosity 
of 10^/6^^ to interpret the results. 

The invariant mass spectrum of the Z pair is shown in fig. 5 for different 
values of the Higgs mass. The background refers to mn oo. The trans- 
verse momentum distribution of the Z bosons is shown in fig. 6. One notices 
that the way the radiative corrections influence the shape of the Higgs res- 
onance differs from the effect shown in fig. 3 for the process tt ^ H ^ zz. 
In particular, the cross section on top of the resonance is slightly increased. 
This is due to interference effects with the nonresonant diagrams in fig. 4. 

4 Conclusions 

We studied the effects of the radiative corrections of enhanced electroweak 
strength on Higgs production by gluon fusion at hadron colliders. The full 
correction to the Higgs resonance shape at next-to-next-to-leading order 
was derived. This was then incorporated in a Monte-Carlo simulation of 
the processes gg — > ZZ and qq ZZ at LHC. 

Compared to the similar radiative corrections to the H ^ VlVl decay, 
the Higgs resonance shape corrections turn out to be relatively small. This 
is due to cancellations between different contributions. Other than for the 
ti ^ H ^ ZZ scattering, in the gluon fusion process gg ^ H ^ ZZ the 
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radiative corrections result in an increase of 10 — 20% of the Higgs signal 
due to interference effects with the nonresonant diagrams. 

We did not address the question of the other Higgs production mech- 
anism, the vector boson fusion. The cross section of this process becomes 
more important as the mass of the Higgs boson increases. For the vec- 
tor boson scattering process, one-loop results are available in the large 
niH limit. Results at two-loop order in the large momentum limit also exist 
p4[ . A complete two-loop analysis would be difficult to perform even via the 
equivalence theorem because one needs to evaluate two-loop box diagrams 
at finite external momenta. This problem deserves further investigation. 
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Figure captions 



Fig.l The radiative corrections of eniianced electroweak strength to the 
process ti ^ H ^ zz. 

Fig. 2 Cut diagrams which contribute to the imaginary part of the Higgs 
selfenergy at three-loop order. These contributions are not included in the 
two-loop H — ^ zz, w'^w~ decay width. 

Fig. 3 Radiative corrections to the Higgs resonance shape in the process 
ti ^ H ^ zz for niH = 850 GeV. The solid line is the squared absolute 
value of the leading order (eq. 2). The thin line corresponds to the NNLO 
expression (eq. 3). 

Fig. 4 Feynman diagrams for pair production by gluon fusion. 

Fj^. 5 Invariant mass distribution of the pairs at LHC. The processes 

considered are gg ZZ — > 2{fi~^iJ,~) and qq ZZ — > 2{ji~^iJ,~). We take 
^/s = 14.5 TeV, and for the outgoing muons we request pT > 20 GeV and 
\yi\ < 2.5. The solid line is the NNLO cross section, the dashed line is the 
tree level cross section, and the dotted line is the background (no Higgs 
production diagram), a) shows the total cross section, and b) shows the 
Higgs signal, with the background subtracted. 

Fj^. ^ Transverse momentum distribution of the Z*^ bosons at LHC. The 
processes considered are gg ZZ 2{p,'^fi~) and qq — > ZZ — > 2(/x"'"/i~). 
We take ^/s = 14.5 TeV, and for the outgoing muons we request pT > 20 
GeV and \yi\ < 2.5. The solid line is the NNLO cross section, the dashed 
line is the tree level cross section, and the dotted line is the background (no 
Higgs production diagram), a) shows the total cross section, and b) shows 
the Higgs signal, with the background subtracted. 
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